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LASER  PROPAGATION  CODE  STUDY 
by 
Edward  B.  Rockower 

I.  INTRODUCTION  and  SUMMARY 

During  the  course  o-f  this  study  a  number  of  laser  propagation 
codes  have  been  assessed  as  to  their  suitability  t or    modeling 
Army  High  Energy  Laser  weapons  used  in  an  ant i -sensor  mode. 
Because  the  Army  battle-field  scenario  requires  nonlinear  1  ?.=  it 
beam  propagation  through  turbulent  atmospheric  conditions, 
including  smoke  and  dust,  these  aspects  have  required  '  ana 
received)  the  greatest  modeling  effort  (see,  e.q.  I i J > 
However,  other  important  features  of  High  Energy  Laser  weapc 
are  the  rectangular  aperture  arid  large  (possibly  up  to  SL)7.) 
central  obscuration  necessitated  bv  the  unstable  resonator 
designs  common  to  most  high  enerqy  lasers. 

The  adequacy  of  the  codes  in  modeling  these,  and  other, 
a  s  p  e  c  t  s  o  -f  o  p  e  r  a  t  i  o  n  a  1  1  a  s  e  r  s  has  r  e  q  u  i  i '  e  d  validation. 
report  will  identify  some  apparent  problems  with  current 
modeling  in  the  codes,  propose  an  interim  "-fix",  where   o  =  sir.'l 
and  outline  an  approach  to  recommenced  further  wor k ,  based  on 
our  study  o-f  the  i  ssues  l  dent  l  f  i  ed  . 

Two  major  categories  of  code  have  been  investigated  in  our 
study: 


1.  Fundamental  wave  optics  codes;  these  codes  start 
from  the  basic  microscopic  laws  of  E-M  radiation  and  implement 
various  phenomenol ogi cal  models  of  the  atmospheric  turbulence 
and  particulates.   A  number  of  schemes  for  numerically 
integrating  the  resulting  approximations  to  the  coupled  partial 
differential  equations  have  been  mechanized  on  computer  'see 
Table  I  for  a  listing  of  these  and  other  codes). 

2.  Moderate  accuracy  system-Level  scaling  law 
codes;  these  codes  start  from  phenomenol ogi cal  and  analytical 
approximations  to  the  E-M  wave  equations,  and  attempt  to 
represent  the  details  of  propagation  through  the  atmosphere  with 
a  few  parameters  (e.g.  a  phase  integral)  based  on  integrated 
properties  of  the  atmospheric  conditions,  light  intensity,  wave 
length,  etc.   These  codes  ars    based  on  experimentally  and 
numerically  (,  +  rom  the  fundamental  wave  optics  codes;  derived 
data  bases. 

A  large  portion  o+  the  analysis  reported  here  is  relative  to 
two  baseline  codes,  one  from  each  of  the  above  categories.    Ln 
the  first  category  is  the  4—D  code  developed  by  Joe  Fleck  and 
Jim  Morris  at  Lawrence  Livermore  National  Laboratory 1 2  J .   in  the 
second,  scaling  law,  category  is  the  bRLPRO  code  developed  by 
Harold  Breaux  of  the  Ballistic  Research  Laboratory,  a  version  ot 
which  is  contained  in  the  HELAWS  code  which  models  Army  high 
energy  ant i -sensor  laser  weapons C 3  J.   During  the  course  of  this 


study,  this  code  has  been  modified  slightly  while  installing  it 
on  the  Naval  Postgraduate  School   (NPS)  IBM  3033  computer  system 
in  order  to  make  it  compatible  with  the  CMS  operating  system  and 
IBM  Fortran  77. 

A  relatively  large  number  of  computer  codes  have  been 
developed  at  various  institutions  in  order  to  model  hiqh  energy 
laser  propagation.   Unfortunately,  the  capabilities  and 
limitations  at  these  cades  a.r<^    not  always  immediatel/  apparent 
from  a  perusal  of  their  document at  1  on .   Some  at  the  limitations 
a.r&    common  to  all  of  the  cades;  e.q.  ,  the  cade  may  onlv  handle 
circular  laser  apertures.  This  limitation  has  not.  stopped  the 
application  o  t-  the  cades  outside  their  range  o'  validity.   Ut 
limitations  are  those,  such  as  assuming  a  vacuum  in  whi.cn  the 
laser  propagates,  which  clearly  preclude  use  of  such  codes  o  h  it 
than  in  a  Space  Wart  are    scenario.    In  order  to  attempt  to  make  a 
preliminary  assessment  of    what  s  available  we  have  relied 
heavily  on  telephone  conversations  with  engineers  and  ana.J  at  = 
at  various  institutions  regarding  their  codes.    In  addii  : 
significant  source  of  information  has  been  a  survey  ot 
e>:  i  --:  t  i  ng  wave-op  tics,  scaling  1  aw  ,  and  si  mp  1  i  f  i  ad  qeometr  i  c 
Laser  propagation  codes  '.see  fable  1)  carried  out  by  Dr.  J  ames 
P.  Reilly  and  co— workers  at  W.  J.  Schafer  Associates,  Inc.  in 
1979  L4I. 


TABLE  1:   PROPAGATION  CODES 

CODE 

TYPE* 

DEVELOPED  BY 

EXERCISED  BY 

APM 

SC 

HEL  Systems  Project  Office 

Army 

BREAUX 

SC 

Ballistics  Research  Laboratory 

Army 

COMBO 

SC 

Air  Force  Weapons  Laboratory 

Air  Force 

EAPM 

SC 

Charles  Stark  Draper  Laboraory 

Charles  Stark  Draper  Lab. 

ESP  III 

SC 

United  Technologies  Research  Center 

Air  Force 

ESP  IIIA 

SC 

United  Technologies  Research  Center 

Air  Force/Army 

ESP  IV 

SC 

United  Technologies  Research  Center 

Air  Force 

GEBHARDT 

SC 

Science  Applications,  Inc. 

Army 

GUTS 

SC 

Air  Force  Weapons  Laboratory 

Air  Force 

HELP  (PROPMD) 

wo 

Air  Force  Weapons  Laboratory 

Air  Force 

JW/4-D 

wo 

Far  Field,  Incorporated 

Far  Field,  Incorporated 

LASE 

SL 

Science  Applications,  Inc. 

Navy 

LASNEX 

GO 

Lawrence  Livermore  Laboratory 

Lawrence  Livermore  Lab. 

LL/SL 

SL 

Lincoln  Laboratory 

Lincoln  Laboratory 

LL/WO 

WO 

Lincoln  Laboratory 

Lincoln  Laboratory 

LLL/4D 

WO 

Lawrence  Livermore  Laboratory 

Lawrence  Livermore  Lab. 

MPLAW 

SL 

Naval  Research  Laboratory 

Navy- 

NOLEC  " 

SC 

Naval  Ordnance  Laboratory 

Navy 

NRL/CHM 

WO 

Naval  Research  Lab/SAI 

Navy 

NRL/JL 

SC 

Naval  Research  Laboratory 

Navy 

PHILLIPS/SL 

SL 

Science  Applications,  Inc. 

Navy 

PSM 

WO 

Charles  Stark  Draper  Laboratory 

CSDL 

SAICOM 

SC 

Science  Applications,  Inc. 

Navy 

SSPARAMA 

wo 

Naval  Research  Laboratory 

Navy 

UTRC/WO 

wo 

United  Technologies  Research  Center 

United  Technologies  Resear ' 
Center 

ZAPM 

SC 

W.  J.  Schafer  Associates,  Inc. 

W.  J.  Schafer  Associates, 

*  WO  -  w. 
SL  -  S 

ave  Optics 
paling  Law 

:ode                   SC  -  Simplify 

GO   -  Geometric 

4 

id  Code 

:  Optics/Hydrodynamic  Code 

II.  APERTURE  SHAPE  PROBLEMS 

In  conversations  with  analysts  and  programmers  regarding 
their  codes  and  in  study  of  the  documention  on  other  codes, 
where  available,  it  is  apparent  that  virtually  all  non— wave- 
optics  codes  force  the  user  to  model  his  laser  as  having  a 
circular  beam  profile  at  the  laser  aperture.   As  mentioned 
above,  most  high  energy  chemical  lasers  have  rectangular 
ap.ertures.   In  addition,  there  is  o+ten  a  central  or  non-centi 
obscuration  of  the  beam  profile  at  the  laser  aperture,  caused  c 
one  mirror  of  the  unstable  resonator  configuration-    in 
attempting  to  exercise  such  a  propagation  code  to  simulate 
rectangular  aperture  lasers,  it  is  reasonable  to  try  to  man 
both  the  same  output  power  as  well  as  output  beam  intensity  .is 
+or    the  real  laser.    If  the  former  condition  is  met,  then  the 
latter  is  equivalent  to  requiring  that  the  anas,    ot-  the  circui 
aperture  (with  or    without  a  central  obscuration)  be  the  same  ?^ 
f-or  the  real  Laser. 

The  following  analysis  is  an  attempt  to  estimate  the 
seriousness  of  the  limitation  to  circular  apertures  and  lead ; 
to  a    proposal  for  an  interim  fix.   Finally,  we  suggest  a 
possible  course  of  further  work  on  this  problem. 

rin  asymptotic  approximation  formula  for  estimating  the 
fraction  of  encircled  energy  within  a  given  radius  has  been 


derived  for  imaging  systems  with  oddly  shaped  apertures C 5] . 
Modulation  Transfer  Function  (MTF)  techniques  were  applied  to  a 
uniformly  illuminated  aperture  having  arbitrary  shape  and 
obscurations,  resulting  in  the  following  formula,  valid  for 
asymptotically  large  values  of  r. 


E(r)  =  1  -  A  f  R/(27T2r) 


where      >         =  laser  wavelength 

f  =  effective  focal  lenqth 

r  =  radial  dimension  in  the  focal  plane 

K'  =  HEL  aperture  per  1  meter  — to— area  ration 

E  =  normalized  encircled  energy,  i.e.  the 

fraction  of  the  energy  transmitted  by  the  aperture  that  falls 
within  a  circle  of  radius  r  about  the  geometrical  focal  point. 

The  main  feature  of  this  result  for  our  purposes  is  that, 
apart  from  the  laser  wavelength  and  system  focal  distance,  the 
traction  of  encircled  energy  depends  only  on  the  ratio  R/r , 
where  R  is  the  ratio  (per 1  meter /area )  for  the  imaging  system 
aperture,  and  r  is  the  radial  dimension  in  the  focal  plane 
within  which  one  wishes  to  determine  the  fraction  at  total 
energy.   Asymptotically,  the  fraction  of  encircled  energy  is 
independent  of  details  of  the  shape,  apart  from  the  value  of  R. 
It  1=  assumed  that  the  laser  beam  propagates  in  a  linear  medium, 


with  no  effects  o-f  atmospheric  f  1  uct»uat  i  ons ,  wind,  etc.  and  that 
the  laser  aperture  is  uniformly  illuminated. 

Analysing  this  result,  we  conclude  that,  at  least  for  the 
restricted  conditions  under  which  it  was  derived,  the  radius  (r  > 
encircling  a  given  fraction  of  energy  increases  linearly  with  R. 
In  other  words,  a  beam  from  an  aperture  with  double  the  /alue  of 
R  will,  asymptotically,  spread  twice  as  much,  from  diffraction. 
Hence,  again  subject  to  the  limitations  of  applicability,  the 
spread  of  a  beam  in  free  space  propagation  will  be  propot-  i  i  i 
to  the  ratio  of  the  lasers  perimeter  to  its  ^res.. 

There  is  a  somewhat  similar  effect  that  is  a  wei J  kn 
consequence  of  imperfections  of  laser  wavefront  qu  ?.lit>  3  '<     tne 
point  where  the  beam  leaves  the  laser  aperture.   Whatever 
cause  of  the  degraded  beam  quality  (phase  front  distar  tions; , 
whether  from  1 nhomoqenei t i es  in  tne  lasinq  medium,  mirror  o; 
lens  imperfections,  etc. ,  the  effects  on  beam  propagation  *\~^ 
represented  with  a  parameter  known  as  the  "beam  quality"  (M). 
The  value  of  M  is  always  greater  than  or"  equal  to  I  and  i_  a 
as  "times  diffraction  limited";  i.e.  ,  the  beam  spread  ii  h  tirnes 
the  diffraction  limited  rate  of  beam  divergence  (proportional  to 
.lambda/D>.   Comparing  the  two  results  for  beam  spread,  we 
identify  a  possible  method  of  compensating,  at  least  partial  . , 
for  rectangular  apertures  and  various  sizes  c-md  types  (e.q.  non- 
central;  of  obscuration.  Our  conclusion  is,  far     rree  space 


propaqation,  that  the  asymptotic  spread  of  lasers  emanating  from 
two  apertures  having  the  same  value  o-f  the  product  li*R  will  be 
the  same.   Hence,  a  better  approximation  to  a  rectangular 
aperture,  possibly  with  a  central  obscuration,  than  just  finding 
the  circular  aperture  with  the  same  total  Are^,     is  to  also 
change  the  value  o-f  beam  quality,  N,  so  that  the  relation 

MR  =  M'R' 

is  preserved. 

For  example,  we  can  model  a  circular  aperture  with  a 
central  obscuration  with  another  circular  aperture  o-f  the  same 
&r<BA,     output  power,  and  with  beam  quality  given  by  M   =  li(R/R'). 
These  two  requirements  lead  to. 
Circular  apertures: 


D'  =  D  /(I  -  F) 


where  F  is  the  fraction  obscuration. 

We  have  also  derived  similar  equations  tor  rectangular  apertures 
with  and  without  obscurations.   The  results  -for  the  latter  ^r& 
presented  in  graphical  form  in  Figure  1.   Each  curve  in  that 
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Figure  1.  Beam  quality  as  a  function  o-f  Length/Width  ratio  (L/W) 
and  Fraction  obscuration  (F). 


figure  gives  the  beam  quality  of  an  "equivalent"  uniform 
circular  aperture  with  the  same  ares,    as  the  rectangular  aperture 
whose  length  to  width  ratio  (L/W)  is  shown  to  the  right  of  the 
curve.   The  curve  for  a  circular  aperture  with  obscuration  F  is 
also  shown. 

It  remains  to  determine  whether  the  above  results  have  any 
relevance  or  usefulness  for  nonlinear  laser  propagation  within  a 
turbulent  atmosphere. 

In  order  to  estimate  the  magnitude  of  the  problem  and  to 
determine  whether  the  similarity  transformation  derived  above 
may  be  useful  we  have  exercised  the  HELAWS  code  (containing 
BRLPRQ)  with  the  following  specific  questions  in  mind: 

t)  How  adequate  is  the  remedy  of  simply  representing  a 
rectangular  aperture  by  using  a  circular  one  of  equal  area? 

2)  How  adequate  is  it  to  represent  a  centrally  obscured 
aperture  by  means  of  a  circui  <Ar  aperture  with  the  same  overall 
3r&a.i     power,  and  (there+ore)  intensity? 

3)  Given  that  we  would  like  to  try  to  represent  different 
aperture  shapes  and  obscurations  with  our  similarity 
transformation,  how  good  a  fix  can  be  obtained  in  the  presence 

°  +  1 i near  atmospheric  effects  such  as  turbulence  and  wind? 
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4)  How  good  is  our  transformation  in  -face  ot  nonlinear 
effects,  i.e.  thermal  blooming? 

5)  How  can  the  existing  codes  be  improved,  possibi     1 
an  extension  ot  the  similarity  trans-formation,  or  otherwise. 

With  re-ference  to  the  above  questions,  we  now  1  ook  at  the 
following  3  figures.   The  data  were  generated  with  the  output 
from  HELAWB ,  using  variations  of  the  base  case  parameters  snown 
in  Table  2. 

On  each  of  the  figures,  curve  <1)     represents  the  results  ror  - 
reference  uniformly  illuminated  circular  aperture  of  diameter 
one  meter,  beam  quality  equal  to  one.   Curve  (2)  presents  the 
results  for  a  uniformly  illuminated  circular  aperture  with  ■:■>. 
central  obscuration  but  with  the  same  a.re'^    and  beam  qua! 
(i.e.  M  -  1)  as  the  base  case.   Since  F  —  .  1,  one  or  our 
previous  formulae  yields  D   -::  1.054*0.   Finally,  curve 
present  5  the  results  for  a  uniformly  illuminated  circular 
aperture  mo  obscuration) ,  equal  area  as  be+ore,  but  with  beam 
quality  calculated  from  our  formula  to  give  the  same  +r<^o     space 
asymptotic  beam  spread  as  the  aperture  with  lcj".  obscur  at  l  on  ■;  1  .  o. 
curve  2)  ,        Our  formula  presented  above  yields  M   =  1  .  38  . 
can  also,  approximately,  be  read  off  ot  figure  4  which  shows  :l  e 
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TABLE  2:  HELAWS  CODE  INPUT  PARAMETERS 

LASER  PARAMETERS 

LASER  CODE R1-B3 

OPERATING  MODE RP 

WAVELENGTH  10.6  MICRONS 

BEAM  QUALITY  (TIMES  DIFFRACTION  LIMITED)  1.39  X  DL 

APERTURE  MIRROR  DIAMETER  1.000  M 

TOTAL  SYS.  JITTER  (1  SIGMA) 10.0  URAD 

VARIABLE  FOCAL  RANGE  EQUAL  TO  TGT  RANGE 

AIMPOINT  OFFSET:  X  COORDINATE 0.00  M 

AIMPOINT  OFFSET:  Y  COORDINATE 0.00  M 

STD.  DEV.  OF  AIM  BIAS  ERROR  (1  SIGMA) .5.0  URAD 

ENERGY  PER  PULSE 10.0  KJ 

PULSE  REPETITION  RATE 5.0  P/SEC 

PULSE  DURATION  (OR  WIDTH) 20.0  USEC 

ATMOSPHERIC  PARAMETERS 

SCENARIO  LOCATION EUROPE 

SCENARIO  LONGITUDE  9.0  DEG 

SCENARIO  LATITUDE. 50.0  DEG 

TIME  OF  YEAR OCT 

TIME  OF  DAY 10.00  HRS 

VISIBILITY  RANGE 7.00  KM 

AMBIENT  TEMPERATURE 10.0  DEG  C 

RELATIVE  HUMIDITY 85.0 

ATMOSPHERIC  PRESSURE  985.0  MB 

TURBULENCE  LEVEL  MODERATE 

REFRACTIVE  INDEX  STRUCTURE  CONST  (AT  1M)  8.40E-14  M**(-2/3) 

SCALING  OPTIONS  EMPLOYED  ...  YES 

MAGNITUDE  OF  WIND  VELOCITY  AT  REF  HEIGHT 1.0  M/S 

REFERENCE  HEIGHT  FOR  WIND  VEL  ALTITUDE  SCALING  1.0  M 

WIND  DIRECTION  ANGLE  (TO  SOUTH=0  DEG,  TO  WEST=90  DEG) 220.  M 

ENGAGEMENT  PARAMETERS 

INITIAL  TARGET  RANGE  4000.  M 

MAXIMUM  TARGET  RANGE 4500.  M 

RANGE  INCREMENT 100.  M 

ALTITUDE  OF  LASER  APERTURE  2.8  M 

ALTITUDE  OF  TARGET 1.5  M 
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required  beam  quality  of  an  equivalent  circular  aperture  laser 
-for  various  shapes  and  values  of  F. 

First,  let's  look  at  -Figures  2  and  3,  showing  the  unbloomed 
1/e  spot  diameter  and  unbloomed  1/e  peak  fluence,  respectively. 
Comparing  curves  (1)  and  (2)  shows  that  simply  trying  to 
represent  a  107.  obscuration  with  a  circular  aperture  of  equal 
area  would  lead  to  gross  errors.   We  now  appeal  to  the  encircled 
energy  theorem  and  note  that  it  is  easy  to  show  that  there  is  a 
rectangular  aperture  CL/W  appro*.  =1.9)  with  the  same  area  and 
value  of  R  as  the  107.  obscured  circular  aperture.   Hence,  by  the 
encircled  energy  theorem  the  rectangle  will  have  the  same 
asymptotic  beam  spread  as  will  the  obscured  circle.   Clearly,  ij 
Dne  tries  to  represent  such  a  rectangle  with  the  unobscured 

ircular  aperture  (curve  1),  then  the  same  gross  errors  will 
result,  even  without  thermal  blooming  (there's  certainly  no 
reason  to  believe  that  the  nonlinear  thermal  blooming  will  help 

matters  here) . 

Still  looking  at  figures  2  and  3,  now  note  curve  (3)  which 
corresponds  to  an    unobscured  circular  aperture  with  the  same 
area  as  the  other  cases  but  with  beam  quality  M  set  by  the 
requirement  that  M*R  matches  that  of  the  107.  obscuration  case. 
We  see  that  the  match  is  excellent  beyond  about  4300  meters. 
Apparently  this  range  defines  the  beginning  of  the  asymptotic 
region  tor  the  conditions  of  these  computer  runs.   Two  things 


o 


>_ 
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are    worth  noting  here:  1)  the  M*R  =  const,  similarity  relation 
appears  to  work  quite  well,  becoming  better  as  the  range 
increases,  lending  credence  to  our  reasoning  regarding  the 
scaling  of  diffraction  beam  spreading,  whether  arisinq  from  M  or 
from  R,  and  2)  remarkably,  our  result  works  well  even  tor 
propagation  in  the  atmosphere,  with  its  turbulence 
f 1 uctuat i ons (descr i bed  as  moderate  for  this  case)  and  wind  (but 
without  thermal  blooming).   Hence,  the  similarity  relation  hold- 
even  beyond  the  range  of  validity  implied  by  the  conditions  of 
its  derivation. 

C  a  n  s  u  1 1  i  n  q  fig  u  r  e  4  ,  1  n  which  t  h  e  t  h  e  r  m  a  1  3.  1  y  b  loomed  l  .  ■  e  5  p  o  c 
diameter  is  plotted,  we  see  that  adding  the  beam  quality  factor 
to  the  equal  area  circular  aperture  (curve  3)  gets  us  only  a 
little  av3r    half  way  to  the  1(3"/.  obscured  aperture  results 
2).   Hence,  when  the  nonlinear  e+fects  of  thermal  blooming  are 
considered,  the  central  obscuration  is  not  so  easily  compensated 
f  o r  b  y  our  beam  q  u  al it y  s  i  m  i 1  a r i t  y  relation.   Ho w e  v  e r .  t  h e  c  o d e 
a  1  r  e  a  d  y  handles  the  10 "'.  o  b  scuration  c  a  s  e  ,  a  n  d  we  a  re  on]  y 
looking  for  a  means  of  extending  the  range  o-i  existing 
propagation  codes  to  include  larger  obscurations  and  rectangular 
apertures.   Clearly,  using  the  principle  of  constant  M*R  qt  < 
p  a r t  w a  y  t  h e re.   Hence,  a  possible  a  p p r o  a c h  i s  t o  q e n e  i 
extend  the  beam  quality  vs.  R  scaling  relation  lor  similarity 
relation'  to  the  nonlinear  region  by  a)  fitting  cm  /es  to 
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existing  code  output,  and/or  b)  having  identified  the  extent  of 
the  problem  and  a  possible  direction  -for  a  phenomenol oqi cal  fix, 
carry    out  studies  with  existing  wave-optics  codes  on  rectangular 
apertures  and  larger  obscurations  to  obtain  the  applicable 
scaling  laws.   Possibly  some  generalization  at  R*M  =  canst, 
which  also  accounts  -for  beam  power,  e.g.  R*M  =  -f  <.  I )  where  i  is 
the  beam  intensity,  will  prove  adequate. 

III.  AIM  BIAS  AND  DISPERSION 

A  number  of  the  codes  surveyed  allow  at  most  two  options 
regarding  focus  of  the  laser  beam.   Either  the  user  must  assume 
exact  focus  on  the  target,  or    the  focus  is  at  infinity.   Whether 
this  is  justified  or  not  should  be  determined  tor  the  specific, 
realistic  pointing  and  tracking,  and  laser  focusing  systems 
being  planned  for.    It  is  not  difficult  to  show  with  some  very 
simple  calculations  that  a  range  error,  whether  bias  (constant 
from  pulse  to  pulse,1  or  dispersion  can  lead  to  appreciable 
spreading  of  the  laser  energy  at  the  true  target  range.   No 
conclusion  can  be  drawn  without  specific  predictions  for  the 
i-ocusing  system  parameters. 

With  regard  to  bias  of  the  laser  system,  a  straightforward 
analysis  of  the  inputs  and  outputs  from  the  HELAWS  code  indicate 
that  there  will  be  cases  in  which  the  low  frequency  atmospheric 
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turbulence  and  wind-induced  aim  bias,  which,  along  with  the 
postulated  levels  of  laser  system  pointing  bias  will  dominate 
the  dispersion  errors.        For  those  cases  with  high  levels  of  i ow 
-frequency  atmospheric  turbulence  and/or  system  aim  bias,  it  mav 
be  that  the  optimal  -firing  will  require  either  some  type  of 
pattern-f lr ing  of  laser  pulses  or  possibly  an  added  dispersion 
(shot  gun  effect)  to  compensate. 

IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

Based  on  ou.r    survey  of  a  number  of  high  energy  laser 
propagation  codes,  it  would  appear  that  none  ana    highiy  precise 
tools  tor  evaluating  Army  laser  weapons  systems  operating  in 
realistic  environments.    That  this  may  b<=    true  is  not  surpri  -in. 
tor  the  wave  optics  codes  which  were  generally  designed  ho  be 
highly  precise  codes,  valuable  for  baselininq  the  more 
p hen omen ol oqi cal  scaling  codes  intended  for  systems  level 
s t u d 1 e s .   Howe v e r ,  t h e y  w ere  u s u all y  n o t  designed  to  include  t  h  ■ 
range  ot  effects  in  an  efficient  computational  tc^r  m  to  be  us<  o 
1 r  1  s y s t  ems  s t  u d 1 e  s . 

However  ,  consider  the  statement  in  L  _'■  J  vol.  1  p.  L-8,  "Reai  Ml 
beams  may  have  rectangular  cross— sect i ons  with  non-central  or 
non— circular  obscurations  with  greater  tor  less1  than  lid/. 
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..  _  However,  given  the  current  primitive  state 

obscuration nowevc  '  y 

of  HEL  weapon  development  and  the  intended  systems  analysis 
applications  of  HELAWS,  the  BRLPRO  assumptions  -for  the  beam 
characteristics  are  sufficiently  accurate."    This  mav  seriouslv 
underestimate  the  errors  in  approximating  a  rectangular  aperture 
or  one  with  larger  than  107.  obscuration  (or  both)  with  a 
suitably  chosen  circular  aperture,  possibly  with  107. 
obscuration.   We  feel  that  our  discussion  in  section  II  of  this 
report  is  sufficient  grounds  to  question  the  above  quote. 

Hence,  it  may  be  advisable  to  consider  some  of  the 
suggestions  outlined  in  section  II,  above,  regarding 
odifications  of  present  codes  to  more  realistically  handle 
ectangular  apertures  and  larger  obscurations.   Ultimately  there 
st  be  some  validation  either  with  a  wave-optics  code  -such  as 


mc 


(liUi 


ore    4-D    code,     the    NRL    SSPARAMA    code,     etc)     or    with 


the  Li  verm 

experiments  on  existing  high  energy  laser  systems. 
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